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Abstract A Monte Carlo (MC) study is performed to evaluate the surface 
tension 7 of spherical membranes that may be regarded as the models of the 
lipid layers. We use the canonical surface model defined on the self-avoiding 
triangulated lattices. The surface tension 7 is calculated by keeping the total 
surface area A constant during the MC simulations. In the evaluation of 7 , we 
use A instead of the projected area Ap, which is unknown due to the fluctuation 
of the spherical surface without boundary. The pressure difference Ap between 
the inner and the outer sides of the surface is also calculated by maintaining 
the enclosed volume constant. Using Ap and the Laplace formula, we obtain 
the tension, which is considered to be equal to the frame tension r conjugate 
to Ap, and check whether or not 7 is consistent with r. We find reasonable 
consistency between 7 and r in the region of sufficiently large bending rigidity 
K or sufficiently large A/N. It is also found that r becomes constant in the 
limit of A/N —>■ 00 both in the tethered and fluid surfaces. 
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1 Introduction 

Contemporary chemical applications of surface models are closely related to 
the study of the lipid bilayers. Initially the lipid membranes were described 
as a two-dimensional fluids with the protein molecules diffusing in them m- 
But the observed variety of the lipid structures in cells [31I1I5] proved that the 
lipid molecules were much more than just a solvent for proteins. Since then 
the systems consisting of such molecules were extensively studied using both 
experimental and computational methods. It should be mentioned that in spite 
of the early Monte-Carlo approach in m, later the researchers mainly used 
the molecular dynamics (MD) technique [3J|11[SJ[7]. Particular attention was 
paid to the relation between the structural and dynamical properties of the 
membranes and their biochemical functions inside the cell. The MD method 
was applied for simulation of lipid bilayers [iiTieisjiiniiiiiiii]. It gave the first 
glance at the membrane structural features at the molecular level. But at 
the same time it proved to be capable of modeling not more than several 
nanoseconds of evolution for a realistic molecular structure in a reasonable 
computational period; the simulation time range from 10 to 60 nanoseconds 
was described in m as a serious achievement. For this reason in [5] even a 
combination of the MD and MC approaches was proposed. 

In the present paper we pay special attention to the final steady configura¬ 
tions of the investigated membranes, and thus we apply the MC technique. We 
study the mechanical characteristics of membranes, such as the bending mod¬ 
ulus, the surface tension 7 and the pressure difference zip =pin—Pout, where 
Pin(Pout) is the pressure inside (outside) the surface. These physical quanti- 

tially reflect the microscopic interactions of the constituent molecules. 

Experimentally the bending modulus and the spontaneous curvature of a 
fluid membrane were obtained using the so-called flicker spectroscopy tech¬ 
nique, in which numerically generated membrane shapes are fitted to exper¬ 
imental ones m- The same method called the contour detection technique 
was applied to experimentally measured shapes of membranes for extracting 
the bending modulus and the surface tension pO] . 

Therefore, it is interesting to study 7 and Zip by means of stochastic sim¬ 
ulation of mechanical processes on the triangulated surface models, such as 
the Helfrich-Polyakov (HP) model [2Tl[22l[2^ . The frame tension r is defined 
via the macroscopic surface energy, which equals rAp, where Ap is called the 
projected area |24] . This Ap may be regarded as the area contained within 
the boundary T(cR^) fixed in R^, while the real surface area is denoted by 
A. In the case of discrete HP model, A is given by the sum of triangle (or 
microscopic) areas. We use the term ’’frame (surface) tension t ( 7 )” for the 
surface tension conjugate to Ap (H) hence force [l 8 llM] . It is widely accepted 
that on the surfaces spanning F the surface tension 7 is in general different 
from the frame tension r [18| . To the contrary, these 7 and t are expected 
to be identical when the surface is sufficiently smooth m and there is no 
difference between A and Ap. 
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In the experimental measurements, the projected area Ap is also used to 
obtain the frame tension r in a variant of the Laplace formula for the cell 
deformed between the parallel plates [251126] . The pressure difference Ap is 
obtained by measuring both the contact area at the plates and the force F 
on the plates, and then the frame tension t may be found using the Laplace 
formula. Thus, the boundary condition imposed on the cell by the plates fixes 
the shape parameters such as the contact area and the radius, which play the 
role of Ap for the spherical membrane. 

However, Ap is not always well-defined on the surfaces without boundary. 
Therefore, the frame tension cannot be obtained directly on the sphere without 
the fixed boundary. Nevertheless, the pressure difference Ap can be obtained, 
and using the Ap and the Laplace formula, we are able to find the effective 
frame tension 7 |y,Lap- 

The purpose of the present study is to find whether or not 7 , 7 |v,Lap and 
the simulation techniques are well defined for spherical membranes without 
the fixed boundary. For this purpose, in the present paper we report the MC 
data including those in [27] with the detailed information on the calculation 
formula in a self-contained manner. During the MC sweeps (MCS), A is kept 
fixed at the constant value Aq (within the error 0.35% for example). Almost 
the same constraint is imposed on the enclosed volume to obtain the pressure 
difference Ap in the constant volume simulations. From Ap and the Laplace 
formula, we obtain the tension 7 |y,Lap, which can be regarded as the frame 
tension r as mentioned above. One additional point that should be noted is 
as follows: The area/volume constant MC simulation (ACMCS/VCMCS) is 
mathematically considered to be connected with the mapping r : M —>■ 
under the constraint of constant area/volume, where M is a two-dimensional 
surface in the context of Polyakov’s model for strings m- These two mappings 
are in general considered to be very different from each other, because the 
surface morphology with the constant area is different from that with the 
constant volume. Thus, it is interesting to see a relation between the surface 
model and these constrained mappings in the simulations results for 7 and 
7|y.Lap- 

It must also be noted on a relation between the fluid surface model de¬ 
fined on dynamically triangulated surfaces and the liquid phase of the lipid 
membranes. The lipid molecules undergo a first order transition, which is the 
so-called main transition, between the gel and liquid phases. The free diffu¬ 
sion of lipids seen in the liquid phase is greatly reduced in the gel phase. On 
the other hand, the triangulated lattice models can also be divided into two 
groups; one is the fixed connectivity model and the other is the dynamically 
triangulated model. The latter one is called the fluid surface model, because 
the vertices (of triangles) can diffuse freely over the surface due to the dynam¬ 
ical triangulation. On the fixed-connectivity lattices, the vertices can fluctuate 
only locally and hence have no free diffusion. Therefore, the model in this pa¬ 
per, which is a fluid model, corresponds to membranes in the liquid phase at 
least. 
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2 Models 

2.1 Surface tension 7 

The discrete Hamiltonian of the model is defined on a triangulated sphere with 
the vertex position ri{i = l, ■ ■ ■, N) and the triangulation T [H] such that 

5'(r,T;K, Ho, Ap) = aSi + KS2 + Upi^iAo) - ApV + Us, 

S'! = -rJ)^ 5*2 = ^(1 - n, • nj), (1) 

ij ij 

where is the Gaussian bond potential and S2 is the bending energy with 
the bending rigidity ull/ksT]] ks is the Boltzmann constant and T is the 
temperature. We should note that T included in S' as a dynamical variable 
is only for the fluid model, and it is fixed for the fixed-connectivity model. 
The Hamiltonian S is considered as the one for the A^-particle system, and 
therefore S can be viewed as the microscopic Hamiltonian. We should note 
that the coefficient a of Si, which is the microscopic surface tension, is fixed 
to a = 1 hence force. The symbol denotes the unit normal vector of the 
triangle i, and of ^2 denote the sum over all nearest neighbor 

vertices and triangles, respectively. The potential 17Fix(Ho), which fixes the 
surface area H to a constant Aq without using the boundary T, is defined by 

TT 1 _ foo (ll-H/Hol > = l/iVy), 

t/Fix(Ho) - I Q (otherwise), 

where Nt is the total number of triangles Nt = 2N — A. The reason why the 
constraint is imposed on the real area A instead of the projected area Ap is 
because A is well-defined while Ap is not. Note also that the surface is allowed 
to have only in-plane deformation if ca is exactly zero, therefore eA must be 
a small positive number. The energy —ApV with the enclosed volume V is 
well-defined only for the self-avoiding (SA) surfaces. Thus, the energy — ApV 
should be included in the Hamiltonian together with the SA potential Us, 
which is defined by 


Us = Y,Us{A,A'), 

AA' 


Us{A,A') 


00 (triangles A A' intersect), 
0 (otherwise). 


(3) 


where "Y^aa' denotes the sum over all non-nearest neighbor triangles. This SA 
potential is considered as an extension of the one in the Doi-Edwards model 
for polymers [28] , and it is slightly simpler as compared to the one in Ref. [29| . 
The partition function Z is defined by 


N 


Z = ]^driexp[-S'(r,r)] 


r 


(4) 
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where denotes the sum over all possible triangulations. This Yt 
eluded in Z only for the fluid surface model; the fixed-connectivity model is 
defined by Z without Yt' prime in f Il^i dvi denotes that the cen¬ 
ter of mass of the surface is fixed at the origin of R^. We comment on the 
phase space measure in Z. The integration measure 0^1 can be replaced 
by dvi{qi/?)Y, where qi is the coordination number of the vertex i and 

a = 3l2 [nnilSI]- Since qf can also be written as gf = exp(crlogg^), this measure 
effectively turns to be the term —a Yi log Qi in tbe Hamiltonian. Therefore, 
this term is expected to influence the distribution of the coordination number 
on the fluid surface model. Therefore the fluid model in general depends on 
the integration measure which includes the coordination number. However, we 
assume tT = 0 for simplicity in this paper just like in Refs. [sniisg. 

The surface tension 7 is mathematically understood in the context of HP 
model m- Indeed 7 is found from the scale invariance of Z, which is expressed 
by the equation 9o, log ^(ar)|a=i = 0. It is easy to see that S 2 and Us are 
scale independent and «S'i(ar) = Q;^5'i(r), ApV[av) = ApV [r). The depen¬ 
dence of S on the variable T, which is necessary for the fluid surface model, 
is not explicitly written for simplicity. Since Z{AQ]av) = Z([/Fix(^o); or) = 
^(C^Fix(o~^ylo); r) = we see that the corresponding derivative 

da logZ{Ao; ar)\a=i can be written as da-^Ao ^ogZ{a~^Ao; r)i9Q,(a“^Ho)|a=i = 
—2 Hoi9ao log^(^o;r) [16]. Thus, we have 


(3{N - l)a3(^-i)-i^(Ho; or) 

E/n dviSi exp [—S'(or)] 

T i 

+3a3(^-i)+2 E/n dviApV exp [—S(ar)] 

T i 


-2Aoa-^ 


dZ{a ^ylojr) 


Z{Ao;ar) 


= 0 , 

oc—l 


(5) 


where S(ar) = a^Si-|-KS 2 -l-C/Fix(Q! '^Ao) — a^ApV + Us, and Yr included 
only for the fluid model. To evaluate the final term in Eq. ®, we should 
recall that the free energy for the surface with the fixed boundary is given by 
F = — \ogZ = TAp+XHei — ApV [13]. Here, ksT {ks and T are the Boltzmann 
constant and the temperature) is set to be ksT = 1 , the symbol Ap denotes 
the projected area of the boundary, and AiJei is the bending energy. However, 
the surface area Aq is fixed while Ap is unknown in the model of this paper, 
because the surface has no fixed boundary, and therefore Ap is replaced by the 
surface area Aq. Because of such replacement, the frame tension r in this case 
should be changed to the surface tension 7 . Thus, we have 


F = - log Z(Ao, V) = jAq + AHei - ApV. (6) 


Note also that the bending energy term AHei does not contribute to the surface 
tension at least in the term 
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—log Z(Ao, F). Moreover, the term AiJei can be neglected from the 
Hamiltonian if the surface area and the volume become sufficiently large (Hq —> 
oo, V —>- 00 ). In this case the terms 7 H 0 and ApV are apparently dominant 
because AiJei is scale independent. For this reason, we do not go into detail of 
the term Ai?ei- Thus, from Eqs. and m we have 

7U = i2{Si)A,Ap - 3 Ap(V}a,Ap - 3iV)/(2Ho), (7) 

where {Si)a,Ap and {V)a,Ap are obtained by ACMCS. 

Here, we should comment on the units of length and energy. Since the unit 
of aSi is [ksT] and the unit of Si is the length squared, we understand that the 
length unit is given by i/ksT/a, where a is the microscopic surface tension. 
Note also that a can be fixed to a = 1 from the scale invariance of the partition 
function Z in Eq. (|3|). Indeed, j' IliLi drt exp [—5'(r, T',a,K, Aq, Ap)] in Z can 
be reduced to fY\f^^dv'^eyi.p[—S{Y',T]l,n,A'Q,Ap')] up to a multiplicative 
constant by the change of integration variable such that r —>■ r' = ^/av (<t7 r = 
{l/y/a)Y'). This scale change does not influence the constant Hq, and therefore 
we introduce the new constant A'^ = aAo for the scaled Hamiltonian. From this 
new Hq, it is easy to see t/Fix(i"^; Ag) = [/Fix(r; Ag) because of the relation that 
A'/Aq = A/Aq and the definition of t/Fix(r; Ag) in Eq. (Ql). It is also easy to 
see that Ap'V = ApV, where Ap' is defined by Ap' = Ap. Thus, using 
the fact that S 2 and Us are scale independent, we prove S{r,T; a, k, Aq, Ap) = 
S{r',T;l,K,A'f^,Ap'). 


2.2 Pressure difference Ap 

The pressure difference Ap can also be calculated using almost the same 
model and procedure as those described above for 7 ]^. Indeed, the constraint 
UpixiAo) in Eq. ([T|) can be replaced by UpixiVo) to fix the enclosed volume V 
to Vo; a-nd the energy — ApV should be removed from the Hamiltonian. In this 
model, Ap is not an input parameter but it is produced as an output. Thus, 
the Hamiltonian is given by 


S{r, T; K, Ho) = 5i + KS 2 + UfM) + Us, ( 8 ) 


where UpixiVo) is defined by replacing AjA^ with H/Vb in Eq. ([2|) 


C/fixCFo) 


00 (|I - H/Vb| > ey), 
0 (otherwise). 


(9) 


In this expression, ey is defined in such a way that the enclosed volume V 
satisfies |H—Vb| <vo, where vq is the volume of the tetrahedron which consists 
of the regular triangles with the area ao = Aq/Nt', Ag =47r (3Vb/47r)^^^ is the 
area of the sphere of volume Vq. We should note that oq in VCMCS is not 
always identical with the mean triangle area A/Nt, because the area A in 
VCMCS is not fixed to Aq. However, it is expected that Aq^A, because the 
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surface always becomes a sphere in VCMCS although the surface fluctuations 
become large except for k—>- oo at least. 

From the scale invariance of the partition function, we have 


Ap\v = {3N-2{Si)v)/{3Vo). 


( 10 ) 


Indeed, we see that da \ogZ{a ^Vq; r)|a=i can be replaced by —log Z(Vo; r) 
[16]. Thus, we have Eq. cni) by using the same procedure as described in the 
previous subsection. The symbol Ap\v is used for the calculated pressure to 
distinguish it from the input parameter Ap in Eq. (fT|). 

It must be emphasized that Ap\v does not suffer from the problem encoun¬ 
tered in the calculation of the surface tension caused by the difference between 
Ap and A, because the energy increment due to the volume variation AV is 
independent of the way how the surface shape changes. On the contrary, the 
increment of the frame tension energy depends only on the change AAp of the 
projected area Ap. This AAp is not always identical to the change AA of the 
microscopic surface area A of the sphere without the boundary. 

3 Monte Carlo Results 

The canonical Metropolis MC technique is used to update the variables r and 
T, where T is updated only for the fluid surface model by using the standard 
bond flip technique. The acceptance rate i?r for r —>• r' = r-hjr can be controlled 
by tuning the radius cq of a small sphere, inside which (5r is chosen randomly. 

We fix ro to maintain i?rRs^Fix —50%, where i?s a-nd i?Fix are the acceptance 
rates for the interactions t/s: t^Fix(^o) or t/Fix(Vb), respectively. The rate i?Fix 
is almost 100% at sufficiently large k while it is 70% ^ 80% in the limit of 
K—>-0. The fact that i?Fix is relatively high implies that the potentials UpixiAo) 
and 17Fix(Vb) for the fixed area and the enclosed volume are suitable to the 
local update MC technique. The rate i?s also depends on k and i?s — 50% or 
more for k—^ 0 while i?s —100% for k—>- 00 . All simulations are performed on 
the surface with N = 1442. The total number of MCS is approximately 5 x 10® 
including 5 x 10® thermalization MCS. 

3.1 Dependence of 7 on the bending rigidity 

Let us check whether the MC results of qjyi are consistent with those of Ap\v. 

We firstly perform VCMCS by varying k with constant Vq; which is fixed by 
using the triangle areas of the initial sphere so that the mean value oq = Aq /Nt 
becomes ao = 0.14, 0.5 (tethered model) and ao=0.12, 0.5 (fluid model). This 
oq approximately equals oq = 0.216 for the initial sphere of size Si/N = 3/2, 
which holds for the equilibrium surfaces without the potential [/Fix(Ao) under 
Ap = 0. From the MC results Ap\v of the tethered and fluid models we obtain 
7|y.Lap using the Laplace formula 


7k.Lap = iR/2)Ap\v, R = (3Vo/4^)^^^ ■ 


(11) 
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(a) (b) 



Fig. 1 (Color Online) 7 vs. k at (a) aQ{= Aq/Nj') = 0.1A (tethered), (b) ao = 0.12 (fluid), 
(c) ao = 0.5 (tethered) and (d) ao = 0.5 (fluid). The symbols Q *5 ^5 V) denote the 

results obtained by the constant volume (area) simulations. The lattice size is A^ = 1442. 


This 7|y,Lap corresponds to the frame tension r. Here we use the symbol 77,Lap 
for T so that one can see how the frame tension is obtained. These results are 
plotted (O) in Figs. [Ha)”(d). Note that the surface shape remains spherical 
in VCMCS with constant Vq even in the limit of k—^ 0 for both tethered and 
fluid models. 

Next, we perform ACMCS with constant area Aq and Ap, both of which 
are supplied by the output of VCMCS that were performed to obtain the data 
denoted by the symbol (O)- Thus, the input parameters for ACMCS are given 
by 

Aq = {A)v (12) 

and 


(i) Ap = -Ap\v, (ii) Ap = 0, 

(hi) Ap = —2Ap\v, (iv) Ap = Ap\v. (13) 

Under these four different conditions in Eq. (TT^ together with the one in Eq. 

we perform ACMCS and obtain the surface tensions 7 ]^ from Eq. ([7]). 
Thus we predict that the results 7 U(k) satisfy 

(i) 77 = 0 (x), (ii) 77 = -77Lap(*), 

(hi) 7 U = 77,Lap(A), (iv) 7 U = -27|v.Lap(v) 


(14) 
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corresponding to the conditions in Eq. (USD, if the surface remains spherical 
(the data are plotted in Fig.[T|). 

The predictions in Eq. da come from the expectations 

^0 = {V)a.Ap, {Si)v = {Si)a,Ap- (15) 

It is easy to prove the predictions in Eq. m using Eq. da. Firstly, performing 
VCMCS with constant Vq we have Ap\v = {5N — 2(5'i)v)/(3Vb) by Eq. (fTUl) . 
Next, the condition (iii) Ap = —2Ap\v is, for example, used as an input as 
well as the constant area Aq = {A)v for ACMCS. Thus, we have from Eq. dZl) 
that 


7U = {2{Si)a,Ap + 6Ap\v{V)a,Ap - 3N) /{2{A)v) 

= ^ - “) 

= (3IV-2(5i)v)/(2(A)v) 

= {R/2) (3N - 2(5i)y) /(3Eo) = {R/2)Ap\v, (16) 

where the relations in Eq. dTHll are used to obtain the third equality. In the 
forth equality, R is defined by R= {3{V)a,Ap)/{A) v which equals {3Vo)/{A)v 
because ofEq. (TTKI) . Thiscan be identified with i?= ( 3 Vb/ 47 r)^/^ in Eq. da, 
because the surface in VCMCS is expected to be a sphere. Thus, we have (iii) 
in Eq. da- The other formulas in Eq. da are also obtained in the same way. 
The basic assumption for Eq. (fHl) is that the surface must be a sphere in both 
ACMCS and VCMCS. 

These equations in Eq. da are almost trivial because the enclosed volume 
Vb is uniquely determined by its surface area Aq if the surface is a complete 
sphere. Thus we understand that the surface tension is consistent to the 
frame tension 7 |y,Lap when the surface is a smooth sphere (in the limit of 
K —oo). This fact is in good agreement with the result of Ref. m that 
the deviation between r and 7 is proportional to the temperature T (at least 
T— >-0), because k oc l/ksT here in this paper. 

Here we show the meanings of the data in Fig. [T] in more detail. Firstly, 
we note that the effective frame tension 7 |y,Lap (O) obtained by VCMCS and 
plotted in Fig. [T] (a) has 7 |v,Lap < 0 (yjy.Lap > 0) for n <3 (k > 3), because 
~ 7 |v,Lap (minus sign) is plotted in the figures. Because of the same reason, the 
data 7 |a (A) obtained by ACMCS and plotted in Fig. [T] (a) has also 7 ]^ < 0 
( 7 !^! >0) for k<3 (k > 3). We can also see that Ap = —2Ap\v, which is the 
one of the input data corresponding to (iii) in Eq. dT51) for this ACMCS, has 
Ap > 0 {Ap < 0) for K < 3 (k > 3) because of Eq. (fTTll . We see from Fig. [T] (a) 
that — 7 U = — 7 |y,Lap, which is just the expected (iii) 7 U= 7 |y,Lap in Eq. (fTTl) . 
is satisfied except for the region for k—>-0. @ 

Figured] (b) shows that the graph of (iii) 7 ]^ = 7 |y 7 ap(A) is slightly broken 
in the region k ~ 5. The reason of this deviation is that the surface shape 
becomes prolate under relatively large negative pressure Ap=—2Ap\v 

at K ~ 5. In the region of small k, such as k < 3, the pressure difference 
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Fig. 2 (Color Online) The ratio ra = ({A)r//47r)^/^/(3Vb/47r)^/®—l VS. K for (a) the tethered 
and (b) the fluid models. The Gaussian bond potential S\/N vs. n for the check of the second 
equation in Eq. (I15l> for (c) the tethered and (d) the fluid models. The symbol ao = 0.5(iii) 
in (c) denotes that both VCMCS and ACMCS are performed with ao = 0.5, and ACMCS is 
performed under the condition (iii) in Eq. II13II . The lattice size is A'" = 1442. 


-2Ap\v, which is used as an input Ap of ACMCS, is positive, and hence the 
surface becomes a sphere in ACMCS. In the region of sufficiently large n the 
surface also becomes a sphere even though Ap = —2Ap\v is negative. Only for 
the intermediate values of k the surface happens to be the prolate because of 
the negative input of Ap. 

We also see a break in the prediction (ii) of Eq. m for the region of 
small K in Figs. [Hc),(d). It comes only from the deviation of Aq from Ap, 
which is unknown, but expected in that region. Indeed, 7 |a (tIasp) is slightly 
smaller (larger) than —7|y,Lap, where 7 |asp is calculated by replacing Aq with 
Asp = 47r (3(V)A,Ap/47r)^^^ in Eq. ([7]). The value of Ap is expected to satisfy 
Asp < Ap < Aq for the condition (ii). Note that the results Ap\v of VCMCS 
become large negative at ag = 0.5 than those obtained at ag = 0.14 and ag = 
0 .12. The large negative input Ap{= Ap\v <0) makes the surface stomatocyte 
[301132] in the ACMCS for ao = 0.5, and therefore the expectation (iv) 7 |a/2 = 
~7|v,Lap is actually not satisfied in both tethered and fluid surfaces in this 
case. These deviations come from the fact that the predictions in Eqs. (HU) are 
only satisfied for a sphere not only in VCMCS but also in ACMCS. Note also 
that the prediction (i) of Eq. (HI, which is not depicted, is correct also in the 
whole region of k including k = 0 for ao = 0.5. 

To show how the surface of VCMCS deviates from a complete sphere, we 
plot the ratio defined by 



































Surface tension and Laplace pressure in triangulated surface models 


11 


(a) (b) 



tethered fluid 


Fig. 3 (Color Online) Snapshots of the fixed-connectivity and the fluid surfaces obtained 
in VCMCS (upper) and ACMCS (lower) at (a) no = 0-14, /t = 0.5, (b) ao = 0.12, k, = 6. The 
scales of the figures are different from each other. 


Tq = ((^)y/47r)^/^/(3Vb/47r)^/^—1 in Figs.[21a),(b), where ((A)y/47r)^/^ is the 
mean radius obtained from the triangle area and (3Vo/47r)^/^ is the one from 
the enclosed volume. We find that ((A)y/47r)^/^ deviates from (3Vb/47r)^/^ 
only in the region of small k both in the tethered and fluid surfaces. The 
reason of the deviation comes from the large surface fluctuations expected at 

K—>-0. 

The second equation in Eq. (fT5l) can also be checked (Figs. [2lc),(d)). The 
symbol ao = 0.5(iii) in Fig. [SJc) denotes that both VCMCS and ACMCS are 
performed with ao = 0.5, and ACMCS is performed under the condition (iii) in 
Eq. (fOl) . The Gaussian bond potential (S'i)y obtained by VCMCS is in good 
agreement with {Si)a,Ap both in the tethered and the fluid surfaces except at 
K—>-0. This equation in Eq. m holds also for the other three conditions in Eq. 
m almost in the whole region of k just like in Figs.[2lc),(d). Since the surface 
is not a complete sphere at sufficiently small k due to the surface fluctuations, 
the equations in Eq. (USD are not always trivial. This leads us to understand 
that the second of Eq. (fTSll implies that the mapping r : M —>■ under the 
constraint of constant area mathematically shares a common property with the 
one under the constraint of constant volume, as mentioned in the introduction. 
In this sense, the simulations for the surface/frame tension in this paper can 
also serve as a check for this property of the constrained mappings. 

Snapshots in Figs. [2Ka),(b) correspond to the data in Fig. [1] for the (a) 
tethered and (b) fluid surfaces. The surfaces in the upper (lower) row are 
those obtained in VCMCS (ACMCS). In ACMCS, Ap is fixed to (a) Zlp= 1.14 
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Fig. 4 (Color Online) — 7|v,Lap (•) and ^\a,Ap (O) ao(= Aq/Nt) for the tethered and 
fluid models, zip = 0.2,—0.3, 0 are assumed for the constant area simulations. The lattice 
size is A^= 1442. 


and (b) Zip = —0.418, both of which correspond to (iii) in Eq. (TT^ . A small 
deviation in (iii) of Eq. (ITil) seen in Fig. [TJb) comes from the fact that the 
surface shape is not spherical but prolate (see Fig. [3Kb)) as mentioned above. 
The surfaces in Figs.|3Ka) are those obtained from the calculations of 7 7 and 
Ap for the checks in Eq. la, where no inconsistency was found for large k at 
least. 


3.2 Dependence of 7 on the area 

Now we proceed to an additional check for the results of 7 |a and 7 |y,Lap ob¬ 
tained in the ACMCS and VCMCS, respectively, by varying Aq/Nt[= oq) 
with fixed k. In contrast to the first check described above, these two differ¬ 
ent simulations can be performed independently. The parameter oq is exactly 
identical to the mean value of the triangle area A/Nt for ACMCS, while oq 
is only an input parameter to fix Vb but almost identical to A/Nt in VCMCS 
in the large n region at least as shwon in Figs. |3Ka),(b). Thus we obtain Z\p|y, 
from which 7 |y,Lap is calculated by means of the Laplace formula in Eq. (Ill|) . 
Therefore we expect that the simulation results shown in Figs.|4Ka)-(f) satisfy 

- 7lKLap(*) = 7U,Zip(0), (17) 

where 7 7,zip is defined by 

iWap = i2{Si)A,Ap - m) /(2Ao). (18) 

The expectation in Eq. (113 can be obtained as follows: The left hand side 
of Eq. (fTTl) is given by using the Laplace formula in Eq. (fTTll such that 


7|y,Lap = -{R/‘^)Ap\v 
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= -{R/2) (3iV-2(5i)y)/(3Fo) 

= -{m-2{Si)v)/{2A^), (19) 

where the second equality comes from Eq. (113, and the final equality is derived 
from the relations i?/(3Vb) = 1 /Aq and 47ri?^ = Aq. These relations are satisfied 
because VCMCS and ACMCS are respectively performed with the constants 
Vb and Aq, which satisfy Aq =47 r (3Vb/47r)^^^. The final expression in Eq. (fT^ 
is identihed with ^\a,Ap in Eq- dT51) using the assumed equation 

{Si)a,Ap = {Si)v, ( 20 ) 

which is understood as a function of oq. This is also a reasonable relation 
expected for the Gaussian bond potentials obtained in both simulations at 
least for sufficiently large ap. The relation in Eq. (EHl) is weaker than those 
of Eq. (1151) in the sense that the surface shape of ACMCS is not necessarily 
identical to the one of VCMCS: stomatocyte, cup-like and dumbbell, even 
branched-polymer, are allowed in ACMCS for the check of Eq. (1171) . The basic 
assumption for Eq. (II3 is that the surface must be a sphere only in VCMCS. 
We note that the surface shape of the latter simulation (=VCMCS) must be 
spherical because the Laplace formula for a sphere is assumed to yield 7 |y,Lap- 
This requirement is always fulfilled as mentioned above. Note also that the 
constant value of yjv.Lap in the limit of oq —>-00 is independent of whether the 
model is tethered or fluid. 

The value ^\a,Ap is identical to the surface tension in Eq. ([71) if Ap = 0. 
However, it must be emphasized that the relation in Eq. (ED does not always 
imply that the surface tension ^\a.Ap can be identified with the frame tension 
dIv.Lap but it only implies that ^\a,Ap is equal to the minus of qjy.Lap obtained 
under the assumed condition. 

An additional comment on the results in Eig.|3]is as follows: The negative 
dlv.Lap comes from the negative output Ap\v obtained in VCMCS (see Eq. 
ED)- In other words, this negative pressure is understood to be the one caused 
by the negative frame tension at sufficiently large oq- This implies that the 
pressure difference becomes zero if the positive pressure Ap = —Ap\v can be 
externally supplied. Therefore, the results in Fig. |7] also imply that the frame 
tension in the limit of oq —^ 00 becomes independent of oq under a suitable 
external pressure Zip(ocl/i?) for spherical surfaces, and that this property is 
also independent of whether the frame tension is negative or not. 

To show how the surface of VCMCS deviates from a complete sphere, we 
plot the ratio 

ro = ((A)v/ 47 r)^/^/( 3 Vb/ 47 r)^/^—1 in Figs.[nDa),(b). We find that ((A)v/ 47 r)^/^ 
is almost identical to ( 3 Vb/ 47 r)^/^ for large oq region. To the contrary, in the 
limit of oq —>-0 the deviations of ((A)y/ 47 r)^/^ from ( 3 Vb/ 47 r)^/^ are relatively 
large for the region of small n. This is also expected from the results shown in 
Figs.[2i;a),(b). 

The equation (1701) can be checked by the ratio rs^ = |1 — (S'i)v/(5'i)yi,4p|, 
which is plotted in Figs. [nDc),(d). The value of Si/N becomes very large in 
the limit of uq —> 00 , and for this reason the ratio is plotted instead of 
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Fig. 5 (Color Online) The ratio ra = ((A)v/47r)^/^/(3Vb/47r)^/^ — 1 vs. ao for (a) the 
tethered and (b) the fluid models in VCMCS. The ratio rg^ = |1 — (Si)y/(S i)A,Ap\ for (c) 
the tethered and (d) the fluid models. The symbol K; = 0.1(Z\p = 0) in (c) denotes that both 
VCMCS and ACMCS are performed with k = 0.1, and ACMCS is performed under Ap = 0. 
The lattice size is N = 1442. 


Si/N. The symbol K = 0.1{Ap = 0) in Fig.jSfc) denotes that both VCMCS and 
ACMCS are performed with a{ = 0.1, and ACMCS is performed under Ap = 0. 
The data in Fig. [S^c) correspond to those in Figs. l4l[b),(c) for the tethered 
model, while the data in Fig.jSJd) correspond to those in Figs. |llje),(f) for the 
fluid model. We see that —>■ 0 (oq —>■ oo) and therefore the relation in Eq. 

(OUll is satisfied at sufficiently large ao in both models. 

Snapshots in Figs. IUa),(b) correspond to data presented in Fig. 0] (c) and 
Fig. 111(f). In ACMCS (lower), Ap is fixed to Ap = —0A (Fig.lHKa)) and Ap = 0 
(Fig. HKb)). The surfaces in Figs. UKa), (b) are obtained from the simulations 
for the calculations of 7 and Ap for the check in Eq. (ITTll . No deviation was 
found in the data from VCMCS and ACMCS for large oq at least. In Eigs.[6Ka), 
(b) the surfaces obtained in ACMCS (Eig. EKb) lower) differ from the sphere, 
however the relation in Eg. (1201) is satisfied and for this reason 7 is considered 
to be calculated consistently with Ap according to the expectation in Eq. (ITT)) , 
although the configurations of the fluid model obtained in VCMCS at large uq 
are fluid only partly (Fig. HKb) upper). 

Note that the acceptance rate of the bond flip R-p is very low {Rp < 0.05) 
in VCMCS for the fluid model at ao(> 0.75). The reason for this is such a 
large energy change as ASi, which is caused by the bond flip; consequently 
the vertex diffusion is localized. For this reason, the surface is not completely 
fluid but - only partly - fluid for large ag in VCMCS for the fluid model. As a 
consequence the density of vertices becomes non-uniform on the surface, and 
therefore the surfaces at large og may be different from the ordinary fluid ones 
at relatively small ag. Indeed, we can see in the snapshot in Fig. EKb) that 
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(a) 


(b) 



tethered 


fluid 


Fig. 6 (Color Online) Snapshots of the fixed-connectivity and fluid surfaces obtained in 
VCMCS (upper) and ACMCS (lower) (a) ao = l, k = 3, and (b) ao = 4, k = 0.2. The scales 
of the figures are different. 

the density is higher in the upper side of the sphere than in the lower side. 
The lower side of the sphere is relatively smooth and is composed of large 
size triangles. Almost all coordination numbers q are expected to be g = 6 . It 
implies that there is no vertex diffusion in the lower side. Thus, the localization 
of vertex diffusion is expected on the surface in Fig. Elb). 


To avoid the localization of vertices, we can introduce the potential U 2 = 
U 2 i{ao) with 



( 21 ) 


where is the area of triangle i and £2 = 0.5 for example. Due to this po¬ 
tential, every triangle area follows (l/2)ao < < (3/2)ao for £2 = 0.5. As 

a consequence, the total area A is also constrained to obey the same relation 
as tti- Recalling that the surface shape is always spherical in VCMCS, we find 
that the total surface area in VCMCS is limited only by the constant volume 
Vo and hence it is not influenced by t/ 2 . Therefore the formula for Aply in 
Eq.dig is not influenced by U 2 as well. The potential U 2 has also no influence 
on the formula for 7 ^ in Eq. o, because the area A is fixed to Aq by UpixiAo) 
in Eq. m in ACMCS. Thus, we can check whether or not 7 ^ is consistent 
with Ap\v using the configurations in which the vertices are prevented from 
localization by the potential 1 / 2 - Indeed, the results in Fig. 01(1) remain almost 
unchanged or relatively close to those in Figs. EKd) and SDe) if the potential 
U 2 is introduced in both ACMCS and VCMCS. In the ACMCS with U 2 the 
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surface configurations are the branched polymer; they are identical to those 
without U 2 , while in the VCMCS with U 2 the vertices freely diffuse over the 
surface even at relatively large oq and hence the configurations are manifestly 
different from those without C/ 2 - 


4 Summary and conclusions 

In this paper, we present the surface tension 7 and the pressure difference Ap of 
spherical membranes together with the calculation techniques on triangulated 
surfaces for the discretization of the canonical surface model of Helfrich and 
Polyakov. Using the real area A of spherical surface, 7 is calculated by means 
of the area constant MC simulations (ACMCS). Using the Laplace formula, we 
evaluate the effective frame tension ylyLap from zip, which is also calculated 
by means of the volume constant MC simulations (VCMCS). Thus, we have 
shown that the surface tension 7 , conjugate to the real surface area A, is 
consistent with 7 |yLap at sufficiently large k or ao(= Aq/Nt)- Our results 
also show that this consistency holds for those data obtained on non-spherical 
surfaces in ACMCS at sufficiently large uq. Thus, the results shown in this 
paper support that the frame tension can be evaluated by 7 by means of 
ACMCS with the proper input data for the constant area and the pressure 
difference under suitable condition such as sufficiently large k or oq. 
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